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ABSTRACT 

Several  acoustic  Doppler  current  profilers  and  vertical  strings  of  temperature,  conductivity,  and  pressure 
sensors,  deployed  on  and  around  the  East  Flower  Garden  Bank  (EFGB),  were  used  to  examine  surface  wave 
effects  on  high-frequency  flows  over  the  bank  and  to  quantify  spatial  and  temporal  characteristic  of  these 
high-frequency  flows.  The  EFGB,  about  5-km  wide  and  10-km  long,  is  located  about  180-km  southeast  of 
Galveston,  Texas,  and  consists  of  steep  slopes  on  southern  and  eastern  sides  that  rise  from  water  depths  over 
100  m  to  within  20  m  of  the  surface.  Three-dimensional  flows  with  frequencies  ranging  from  0.2  to  2  cycles  per 
hour  (cph)  were  observed  in  the  mixed  layer  when  wind  speed  and  Stokes  drift  at  the  surface  were  large.  These 
motions  were  stronger  over  the  bank  than  outside  the  perimeter.  The  squared  vertical  velocity  w2  was 
strongest  near  the  surface  and  decayed  exponentially  with  depth,  and  the  e-folding  length  of  w2  was  2  times 
larger  than  that  of  Stokes  drift.  The  2-h-averaged  w2  in  the  mixed  layer,  scaled  by  the  squared  friction  velocity, 
was  largest  when  the  turbulent  Langmuir  number  was  less  than  unity  and  the  mixed  layer  was  shallow.  It  is 
suggested  that  Langmuir  circulation  is  responsible  for  the  generation  of  vertical  flows  in  the  mixed  layer,  and 
that  the  increase  in  kinetic  energy  is  due  to  an  enhancement  of  Stokes  drift  by  wave  focusing.  The  lack  of 
agreement  with  open-ocean  Langmuir  scaling  arguments  is  likely  due  to  the  enhanced  kinetic  energy  by  wave 
focusing. 


1.  Introduction 

The  oceanic  surface  layer  is  controlled  by  momentum 
and  buoyancy  fluxes  near  the  surface  and  by  entrain¬ 
ment  fluxes  at  the  base  of  the  mixed  layer.  These  fluxes 
generate  a  wide  range  of  physical  processes  between  the 
air-sea  interface  and  thermocline,  and  play  a  vital  role  in 
climate  variability,  biological  productivity,  and  marine 
pollution.  In  the  highly  unstable  boundary  layer,  buoy¬ 
ancy  dominates  the  production  of  turbulent  kinetic  en¬ 
ergy  (TKE)  (e.g.,  Shay  and  Gregg  1986;  Brubaker  1987; 
Lombardo  and  Gregg  1989).  The  production  of  TKE 
driven  by  surface  cooling  initially  modifies  the  vertical 
velocity  component  and  is  then  transferred  by  turbulent 
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pressure  forces  into  horizontal  velocity  components, 
thus  generating  coherent  flows  in  geophysical  boundary 
layers  (see,  e.g.,  Lumley  and  Panofsky  1964;  Thorpe 
2005).  When  the  surface  buoyancy  flux  is  small  and  the 
boundary  layer  is  near  neutral,  the  production  of  TKE  is 
dominated  by  shear  production  driven  by  surface  wind 
stress  (Dillon  et  al.  1981;  Soloviev  et  al.  1988).  Shear 
production  modifies  horizontal  velocity  fluctuations, 
especially  in  the  downwind  direction  before  transferring 
energy  to  the  crosswind  and  vertical  energy  components 
via  turbulent  pressure  fluctuations.  These  ideas  have 
been  used  in  turbulent  closure  schemes  in  the  oceanic 
mixed  layer  (e.g.,  Mellor  and  Yamada  1982). 

However,  transfer  of  momentum  and  buoyancy  at  the 
air-sea  interface  is  complicated  by  the  presence  of  sur¬ 
face  waves.  Surface  waves  boost  the  vertical  transport  of 
surface  fluxes  through  (i)  enhancement  of  TKE  dissipa¬ 
tion  by  breaking  and  (ii)  generation  of  production-scale 
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coherent  structures  in  the  mixed  layer  via  wave-mean 
flow  interactions.  Several  investigators  (Agrawal  et  al. 
1992;  Anis  and  Mourn  1995;  Drennan  et  al.  1996;  Melville 
1996;  Terray  et  al.  1996;  Soloviev  and  Lukas  2003)  report 
that  the  dissipation  rate  of  TKE  in  the  wave  boundary 
layer  is  up  to  an  order  of  magnitude  larger  than  the  pre¬ 
dictions  from  wall  boundary  layer  scaling.  The  enhanced 
TKE  dissipation  rate  by  wave  breaking  is  an  important 
factor  for  momentum,  heat,  and  gas  transfer  rates  at  the 
air-sea  interface. 

As  a  response  to  the  wind  over  the  ocean,  roll  vortices 
are  formed  roughly  parallel  to  the  wind  direction.  These 
circulation  patterns  can  be  recognized  by  the  lines  of 
convergence  at  the  surface,  as  noted  by  Langmuir 
(1938),  and  they  are  referred  as  “Langmuir  cells.”  The 
orientation  of  these  Langmuir  cells  is  roughly  opposed 
to  the  orientation  of  wind-/shear-driven  turbulent  eddies 
in  the  mixed  layer.  Craik  and  Leibovich  (1976)  describe 
a  plausible  theory  for  the  generation  of  roll  vortices.  The 
Craik-Leibovich  mechanism  requires  interaction  between 
surface  waves  and  wind-driven  shear  flows  to  generate 
a  forcing  function  governed  by  the  wave-induced  Stokes 
drift  and  vertical  shear  of  mean  currents.  In  Langmuir 
turbulence,  the  Stokes  production  term  generates  TKE  in 
both  crosswind  and  vertical  directions.  Detailed  discus¬ 
sions  of  theory  and  relevant  references  are  in  Leibovich 
(1983)  and  Thorpe  (2004). 

Several  investigators  observed  convergence  zones 
on  the  surface  layer  and  the  associated  crosswind  and 
vertical  velocities  in  surface  convergence  zones  (Weller 
and  Price  1988;  Smith  et  al.  1987;  Pollard  and  Thomas 
1989;  Zedel  and  Farmer  1991;  Smith  1992,  1998; 
Plueddemann  et  al.  1996).  Weller  and  Price  (1988)  re¬ 
port  three-dimensional  structure  in  mixed  layer  velocity 
and  narrow  regions  of  downwelling  flow  within  the 
mixed  layer,  in  coincidence  with  bands  of  convergent 
surface  flow.  They  further  report  that  the  downwelling 
flow  in  the  middle  of  the  mixed  layer  was  up  to  0.2  m  s-1, 
and  that  the  downwind,  downwelling  flow  had  a  jetlike 
structure.  Coherent  vortices  in  the  mixed  layer  have  been 
studied  using  Doppler  sonar  techniques  (e.g.,  Zedel  and 
Farmer  1991;  Plueddemann  et  al.  1996;  Smith  1992, 1998) 
and  these  investigators  report  the  consistency  of  wind/ 
wave  forcing  of  Langmuir  circulation.  In  the  existing  lit¬ 
erature  on  Langmuir  circulation  studies,  the  majority  of 
the  observations  is  from  deep  water  environments  and 
there  are  only  a  few  references  related  to  shallow-water 
environments  (e.g.,  Hunter  and  Hill  1980;  Marmorino 
et  al.  2005;  Gargett  and  Wells  2007).  In  a  shallow-water 
environment,  the  physical  setup  is  different  since  top 
and  bottom  boundary  layers  can  merge  under  wind  and 
buoyancy  forcing,  and  thus  bottom  frictional  effects  can 
become  a  significant  factor  in  mixed  layer  dynamics. 


Gargett  and  Wells  (2007)  report  that  Langmuir  circula¬ 
tion  plays  an  important  role  in  sediment  transport  in  the 
shallow-water  environment. 

There  are  numerous  numerical-modeling  studies,  mainly 
large-eddy  simulations  (LESs),  to  examine,  quantify,  and 
scale  Langmuir  turbulence  in  the  oceanic  mixed  layer 
(Skyllingstad  and  Denbo  1995;  McWilliams  et  al.  1997; 
Skyllingstad  et  al.  1999;  Li  et  al.  2005;  Harcourt  and 
D’Asaro,  2008;  Grant  and  Belcher  2009).  The  papers 
cited  here  are  just  a  few  references  out  of  many  in  the 
literature  on  this  subject.  In  the  majority  of  these  models, 
Langmuir  circulation  is  simulated  by  specifying  Craik  and 
Leibovich  (1976)  vortex-force  parameters  while  specify¬ 
ing  a  Stokes  drift  corresponding  to  a  monochromatic 
surface  wave  field.  These  studies  describe  the  formation 
of  Langmuir  circulation,  scaling  of  Langmuir  turbulence, 
and  the  distinction  among  convection,  shear-driven  tur¬ 
bulence,  and  Langmuir  flows  in  the  upper  ocean. 

However,  there  are  still  unresolved  issues  such  as 
quantitative  understanding  of  the  three-dimensional 
structure  of  Langmuir  circulation  and  its  interaction 
with  convective-  and  shear-driven  mixing  and  with  in¬ 
ternal  waves  in  the  thermocline,  and  scaling  laws  for 
Langmuir  turbulence.  The  main  objectives  of  this  study 
are  to  examine  surface  wave  effects  on  high-frequency 
flows  over  the  East  Flower  Garden  Bank  (EFGB),  to 
quantify  spatial  and  temporal  characteristics  of  these 
high-frequency  flows,  and  to  examine  the  scaling  laws 
of  vertical  velocity  based  on  Langmuir  circulation  the¬ 
ory.  In  this  paper,  we  describe  measurements  of  high- 
frequency  currents  over  the  EFGB.  The  bank  is  located 
at  the  edge  of  the  Texas-Louisiana  continental  shelf 
and  rises  from  depths  of  100-150  m  to  about  20-m  below 
the  surface.  When  compared  with  previous  observations 
from  shallow  (e.g.,  Gargett  and  Wells  2007)  and  deep 
water  (e.g.,  Smith  1998),  our  measurements  represent 
both  shallow-  and  deep-water  conditions.  The  vertical 
velocity  fluctuations  described  here  show  strong  corre¬ 
lations  with  winds  and  surface  wave  properties.  We 
suspect  that  these  velocity  fluctuations  are  a  manifesta¬ 
tion  of  Stokes  drift-induced  wave-current  interactions 
such  as  Langmuir  circulations  during  wind  events.  Here, 
we  examine  high-frequency  velocity  fluctuations  in¬ 
cluding  temporal,  lateral,  and  vertical  distributions, 
and  characteristic  vertical  length  scales  for  wide  ranges 
of  wind  speeds  and  Stokes  drifts.  The  paper  is  orga¬ 
nized  as  follows.  Section  2  describes  platforms  and 
measurements.  Section  3  describes  observations  of  winds, 
surface  wave  statistics,  and  background  currents  and 
hydrography.  High-frequency  velocity  fluctuations  and 
scaling  of  vertical  velocity  are  described  in  section  4. 
A  discussion  is  given  in  section  5.  Finally,  major  find¬ 
ings  are  summarized  in  section  6. 


August  2013 


WI JESEKER A  ET  AL. 


1629 


Fig.  1.  Location  of  the  Flower  Garden  Banks  National  Marine  Sanctuary.  Red  circles, 
representing  salt  domes  or  banks,  are  shown  with  selected  bathymetric  contours.  The  black  and 
white  inset  illustrates  bathymetric  contours  of  EGFB  and  the  locations  of  bottom-mounted 
ADCP  (M1-M5;  filled  circles)  and  string  (S1-S4;  filled  stars)  moorings.  The  contour  interval  is 
5  m.  The  filled  triangle  denotes  NDBC  buoy  42047. 


2.  Platforms  and  measurements 

Closely  spaced  observations  of  velocity  (i.e.,  u ,  v ,  w), 
temperature  T,  conductivity  C,  and  pressure  P  were 
collected  between  December  2010  and  December  2011 
over  the  EFGB  as  part  of  projects  sponsored  by  the 
Naval  Research  Laboratory  (NRL)  [Mixing  Over  Rough 
Topography  (MORT)]  and  by  the  Bureau  of  Ocean 
Energy  Management  (BOEM)  [Currents  over  Banks 
(COB)].  The  year-long  observational  program  consists  of 
two  6-month  mooring  deployments  (from  December 
2010  to  June  2011  and  from  June  to  December  2011),  and 
a  2-week  intensive  experiment  from  26  May  to  14  June 
2011.  A  detailed  discussion  of  instrumentation,  data  col¬ 
lection,  sampling  methods,  and  data  processing  is  in 
Teague  et  al.  (2013,  hereafter  TG)  and  therefore  we  just 
describe  a  subset  of  the  data,  relevant  to  this  study,  which 
is  mainly  from  measurements  made  during  the  second 
half  of  the  experiment  (June-December  2011). 

Currents  and  hydrographic  fields  were  collected  from 
five  trawl-resistant  bottom-mounted  ADCP  moorings 
(M1-M5),  referred  to  as  “Barnys”  because  of  their 
barnacle-shaped  dome  (Perkins  et  al.  2000),  and  four 


subsurface  string  moorings  (S1-S4),  equipped  with 
temperature,  conductivity,  and  pressure  (TCP)  sensors. 
Bathymetry  and  locations  of  moorings  at  the  EFGB  are 
shown  in  Fig.  1.  M1-M4  and  S1-S4  were  deployed  at  the 
four  corners  of  the  bank  at  water  depths  of  98, 105, 127, 
and  106  m,  respectively,  and  M5  was  deployed  at  a  water 
of  depth  of  47  m,  about  3-km  north  of  the  bank  peak 
(about  20-m  depth)  (Fig.  1).  The  separations  between 
Barnys  and  associated  string  moorings  were  about 
200  m.  The  moorings  were  equipped  with  Teledyne  RD 
Instruments  Workhorse  ADCPs  operating  at  600  (M5) 
and  300  kHz  (M1-M4)  with  four  transducers  each  with 
an  incidence  angle  of  20°.  ADCPs  recorded  near-full 
water  column  current  profiles  every  15  min  (an  ensem¬ 
ble  average  of  120  realizations).  The  ADCPs  at  Ml,  M2, 
M4,  and  M5  had  2-m  vertical  resolution  and  at  M3  had 
4-m  vertical  resolution.  Shallowest  ADCP  bins  at  Ml, 
M2,  M3,  M4,  and  M5  were  7.5-,  10.4-,  16.4-,  11.5-,  and 
5.5-m  below  the  mean  water  line,  respectively,  while  the 
deepest  bins  were  off  the  bottom  by  6.5,  4.6, 11,  4.5,  and 
3.5  m,  respectively.  The  standard  deviations  for  15-min- 
averaged  currents  were  less  than  0.5%  of  the  water  ve¬ 
locity  at  the  five  mooring  sites.  Because  moorings  are 
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bottom  mounted,  the  platform  motion  is  not  a  factor 
for  velocity  measurements.  Apart  from  water  column 
velocities  at  M5,  surface  wave  parameters  such  as  wave 
height  spectra  and  directional  wave  spectra  were  col¬ 
lected  from  the  ADCP  (Strong  et  al.  2000).  The  spectral 
distribution  of  wave  energy  E(f  \  6)  in  frequency  /  and 
direction  6  was  provided  by  a  bottom-mounted  upward- 
looking  ADCP  (600  kHz)  based  on  the  Teledyne  RD 
Instrument’s  Wave  Array  Technique.  It  treats  the 
ADCP  depth  cell  bins  as  an  array  of  independent  sensors 
for  acquiring  wave  orbital  velocities.  The  frequency- 
direction  spectrum  E(f,  6)  is  then  derived  using  an  array¬ 
processing  algorithm  for  velocities  acquired  every  2  h  for 
5  min  at  a  2-Hz  sampling  rate.  Three  Barnys,  at  M2,  M4, 
and  M5,  were  also  equipped  with  high-resolution  pres¬ 
sure  sensors  (Ppods)  developed  by  Mourn  and  Nash 
(2008).  The  other  two  Barnys  (Ml,  M3)  contained  SeaBird 
Electronics  wave/tide  gauges.  The  sampling  rates  of  the 
Ppods  were  1  Hz  and  the  resolutions  were  about  0.14  mm 
orlPa  (1  Pa  =  lNm~2). 

Four  TCP  string  moorings  were  deployed  at  M1-M4 
during  the  first  half  of  the  deployment  and  three  were 
redeployed  at  M2-M4  during  the  second  half.  Each 
string  mooring  contained  8-12  instruments  that  were 
approximately  equally  spaced  between  7  and  12  m  below 
the  surface  and  1-m  above  the  bottom  (TG).  Some  of 
the  instruments  recorded  only  T  and  C  while  others  re¬ 
corded  T,  C,  and  P.  The  large  majority  of  the  instruments 
were  SeaBird  Electronics  MicroCats  (SBE37).  A  limited 
number  of  Aqua  Trolls  that  recorded  7,  C,  and  P  were 
distributed  on  each  of  the  string  moorings.  The  Micro- 
Cats  sampled  every  6  min  with  accuracies  of  ±0.002°C 
and  3  /jlS  cm-1.  The  Aqua  Trolls  sampled  every  12  min 
with  lower  accuracies  of  ±0.1°C  and  ±0.5%  of  the  data 
magnitude  plus  1/rScm-1.  The  conductivity  sensors 
sometimes  failed  toward  the  ends  of  the  records  as 
a  result  of  biofouling.  We  do  not  address  processing 
and  data  quality  issues  here,  because  those  topics  are 
discussed  in  TG. 


3.  Winds,  waves,  and  background  currents 

a.  Wind  field 

Wind  observations  were  collected  hourly  at  the 
southern  edge  of  the  EFGB  at  the  National  Data  Buoy 
Center  (NDBC)  station  42047  (Fig.  1)  and  at  NDBC 
station  42046  located  approximately  44  km  west  of  the 
EFGB.  Winds  at  42046  were  used  when  42047  was  out  of 
service,  which  was  about  half  of  the  year-long  mooring 
period.  Winds  at  two  buoy  locations  were  similar  when 
they  were  available.  No  wave  data  were  available  at 
these  sites.  For  this  study  wind  measurements  were  from 


42046.  The  buoy  provides  wind  velocity  (U4)  at  4  m  (i.e., 
Z4)  above  the  sea  surface.  The  standard  10-m  (i.e.,  Z10) 
wind  velocity  (U10)  was  estimated  iteratively  by  as¬ 
suming  that  the  near-surface  atmospheric  boundary 
layer  was  near  neutral,  where 

u10  -  u4  +  U**-1  log (Z10/Z4) ,  (la) 

U*|U*|  =T/pw,  (lb) 

T  =  PflO>^  lolUol’  ani^  (lc) 

CD  =  10-4(— 0.016|U10|2  +  0.967|U10|  +  8.0558) 

(Id) 


(Hwang  2011).  The  friction  velocity  is  U*,  k  is  the  von 
Karman  constant  (0.4),  t  is  the  wind  stress  at  the  sea 
surface,  CD  is  the  drag  coefficient,  pw  is  the  density  of 
seawater,  and  pa  is  the  density  of  air.  Winds  ( U\ 0)  for 
September-December  were  dominated  by  north-south 
changes  in  direction  caused  by  passages  of  numerous 
fronts  (Fig.  2a).  Winds  from  May  to  August  were  primarily 
out  of  the  southeast.  Wind  events  with  magnitudes  larger 
than  8 ms-1  lasted  from  1  to  6  days.  Wind  events  as¬ 
sociated  with  passages  of  cold  fronts  were  limited  to 
1-3  days. 

b.  Surface  waves  and  Stokes  drift 

The  wave-induced  Stokes  drift  profile  for  a  random 
wave  field  U$(z)  is  calculated  as  an  ensemble  average  of 
wave  components  (Kenyon  1969),  where 


U5(z)  =  4tt 


f/c  ftlT 

j  jQ  fkE(f,0)G(f,z)dOdf. 


(2) 


The  function  G(/,  z)  is  expressed  as  G(f ,  z)  = 
cosh[2k(H  +  z)][sinh(2kH)]-1,  where  k  =  |k|  is  the 
wavenumber  and  2irf  =  [g/c  tanh(/cH)] 1/2 .  The  water 
depth  is  H ,  the  gravitational  acceleration  (9.81ms-2) 
is  g,  and  the  vertical  depth  coordinate  (negative  down¬ 
ward)  from  the  mean  sea  surface  is  z-  The  ADCP  wave 
measurement  system  provides  E(f,0)  for  a  frequency 
band  varying  from  minimum  frequency  /min  to  cutoff 
frequency  fc.  This  frequency  band  varies  with  environ¬ 
mental  conditions  and  during  the  experiment:  /min  = 
0.0098  Hz  and  fc  =  0.3535  Hz.  Stokes  drift  was  computed 
at  M5  where  H  =  47  m. 

Figure  3  shows  the  ratio  of  surface  Stokes  drift  Us0  = 
|Us(z  =  0)|  to  Uio  as  a  function  of  U10.  In  general,  the 
direction  of  estimated  Stokes  drift  follows  the  direc¬ 
tion  of  the  wind  closely.  Data  presented  in  Fig.  3  are 
within  15°  of  the  wind  direction.  The  scatterplot  shows 
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Fig.  2.  Time  series  between  2  Jun  (year  day  152)  and  8  Dec  2011  (year  day  342).  (a)  Wind 
speed  Uio  (blue)  and  direction  QUvo  (green),  (b)  Surface  wave  height  Hs  (blue)  and  direction  of 
wave  corresponding  to  peak  wave  energy  9P  (green),  (c)  Frictional  velocity  U*  (blue)  and 
Stokes  drift  velocity  at  surface  UsQ  (green),  (d)  High-pass-filtered  pressure  (Row)  at  47-m  depth 
(bottom)  from  Ppod  pressure  sensor  at  M5;  Pgw/1000  is  plotted,  (e)  Time  series  of  15-min- 
averaged  vertical  velocity  at  5.5-m  depth  at  M5,  where  water  depth  is  47  m.  (f)  Depth-averaged 
squared  vertical  velocity  (w2H)  at  M5.  Targeted  wind  events  are  numbered  1-12. 


2-h  estimates  of  the  fractional  Stokes  drift  Us0/U±o  and 
the  bin-averaged  Us0/U\q  over  1ms-1  wind  speed  in¬ 
tervals  with  95%  confidence  limits  (circles  with  bars).  As 
a  reference,  we  also  show  Us0/U\o,  derived  from  di¬ 
rectional  wave  spectra  and  wind  speeds  from  a  nearby 
NDBC  directional  wave  buoy  (42019)  located  170-km 
east  of  the  EFGB  (Fig.  3).  Earle  et  al.  (1999)  showed 
that  the  directional  wave  distribution  can  be  derived 
from:  the  directional  Fourier  series  method,  the  Maxi¬ 
mum  Likelihood  Method  (MLM),  and  Maximum  En¬ 
tropy  Method  (MEM).  In  this  study,  the  MLM  was  used 
to  compute  surface  Stokes  drift  from  the  NDBC  buoy 
data.  We  adapted  this  technique  to  be  consistent  with 
the  Iterative  Maximum  Likelihood  Method  used  for 
deriving  AD  CP  directional  wave  spectra  (Strong  et  al. 
2000).  We  also  computed  surface  Stokes  drift  from  the 
Fourier  method,  and  these  results  are  plotted  in  Fig.  3. 
The  fractional  Stokes  drift  at  M5  increases  from  about 
0.1%  to  about  0.6%  as  wind  speed  increases  from  5  to 


15 ms-1.  For  winds  greater  than  10ms_1,  the  ratio 
USo/U10  at  M5  is  about  0.6,  while  that  at  the  nearby 
NDBC  buoy  is  about  0.7  for  the  MLM,  and  is  about  0.8 
for  the  Fourier  method  (Fig.  3).  Note  that  the  cutoff 
frequency  of  NDBC  observation  is  0.485  Hz. 

Both  wave  measurement  systems  and  wave-simulation 
models  impose  cutoff  frequencies  because  of  obser¬ 
vational  constraints  and  computational  limitations. 
Therefore,  wave  spectral  components  below  this  cutoff 
frequency  (i.e.,  fc)  are  available  for  the  estimate  of 
Stokes  drift  [Eq.  (2)].  Contributions  from  high-frequency 
spectral  components  can  be  approximated  by  assuming 
a  spectral  tail  shape  and  directional  spreading  function 
for  frequencies  higher  than  fc  (Webb  and  Fox-Kemper 
2011;  Rascle  et  al.  2008;  Belcher  et  al.  2012).  The  spectral 
tail  with/-5  contributes  30%  to  the  surface  Stokes  drift 
compared  to  the  truncated  spectrum  with  fc  =  0.41  Hz 
(Belcher  et  al.  2012).  Rascle  et  al.  (2008)  show  that 
the  contribution  can  be  about  0.4%  of  wind  speed  for 


1632 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Volume  43 


Fig.  3.  Ratio  of  surface  Stokes  drift  to  wind  speed,  UsJUio  vs  U\o.  Red  circles  are  bin- 
averaged  UsJUio  for  1  ms-1  wind  speed  intervals.  The  error  bars  show  one  std  dev.  The  ma¬ 
genta  and  green  dashed  lines  represent  bin-averaged  results  from  NDBC  buoy  42019  for  all 
of  2011  based  on  directional  spectra  derived  from  the  MLM  and  Fourier  series  methods, 
respectively.  The  black  dashed  line  represents  a  value  of  1.62%  derived  from  the  PM  spectral 
model. 


simulated  wave  spectra  with  f~5  slope  and  a  constant 
directional  spreading  at  frequencies  beyond  fc  =  0.41  Hz. 
The  computation  of  the  high-frequency  contribution  re¬ 
mains  uncertain  because  of  the  lack  of  information  about 
the  high-frequency  directional  wave  spectrum. 

It  is  noted  that  Stokes  drift  resulting  from  shorter 
waves  penetrate  into  a  very  shallow  depth  near  the  ocean 
surface.  While  short-wave  breaking  in  a  shallow  layer  is 
important  for  mixing  near  the  surface  (e.g.,  McWilliams 
et  al.  1997),  the  contribution  of  high-frequency  (short) 
wave  forcing  to  Langmuir  turbulence  is  not  clear 
(Belcher  et  al.  2012).  Therefore,  in  the  following  analysis, 
we  limit  our  computation  of  Stokes  drift  to  the  imposed 
cutoff  frequency  similar  to  Kukulka  et  al.  (2011).  A  de¬ 
tailed  discussion  of  cutoff  frequency  effects  on  the  Stokes 
drift  estimate  is  beyond  the  scope  of  this  study. 

As  illustrated  in  Fig.  3,  the  ratio  UsJUio  is  a  constant 
(=  0.0162)  for  the  Pierson  and  Moskowitz  (1964,  here¬ 
after  PM)  spectrum  (Li  and  Garrett  1993;  Webb  and 
Fox-Kemper  2011).  The  PM  spectrum  represents  a  fully 
developed  and  unidirectional  wave  field.  However,  the 
ratio  UsJUio  is  shown  to  be  small  relative  to  the  PM 
spectral  model  and  to  vary  with  sea  state  development  for 
spectral  models  of  developing  seas  (Webb  and  Fox- 
Kemper  2011;  Harcourt  and  D’Asaro  2008). 


Time  series  of  significant  wave  height  (i.e.,  Hs ),  di¬ 
rection  (i.e.,  Op)  of  dominant  waves,  and  magnitude  of 
Stokes  drift  at  the  surface  (i.e.,  Us0),  derived  from 
ADCP  wave  measurements  are  shown  in  Figs.  2b  and  c. 
During  wind  events,  wave  height  was  about  3  m  over  the 
bank  and  the  dominant  wave  period  was  about  8  s. 
Stokes  drift  was  about  4-8  cm  s_1,  which  was  a  factor  of 
2-4  times  larger  than  the  friction  velocity  induced  by 
surface  winds  (Figs.  2a-c).  In  most  cases,  the  dominant 
surface  waves  propagated  in  the  direction  of  the  wind. 
Surface  wave-produced  pressure  at  the  bottom  (Fig.  2d) 
was  measured  by  the  Ppod  in  the  M5  Barny  mooring 
(47-m  depth).  Surface  waves  with  a  dominant  wave  pe¬ 
riod  of  8  s,  from  bottom-pressure  measurements,  are 
consistent  with  ADCP  wave  measurements  (Fig.  2b). 
The  timing  and  the  duration  of  the  wave  signal,  mea¬ 
sured  from  both  bottom-pressure  and  ADCP  surface 
wave  measurements,  agreed  remarkably  well. 

c.  Mixed  layer  depth  and  background  currents 

The  stratification  was  strong  and  mixed  layer  depth 
(hereafter  referred  as  MLD  or  subscript  D)  was  shallow 
(—10-15  m)  during  the  summer  months  (Fig.  4a).  The 
stratification  became  weak  and  the  mixed  layer  deep¬ 
ened  during  the  passage  of  atmospheric  fronts  in  fall  of 
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Fig.  4.  (a)  The  color  shading  shows  the  time-depth  temperature  (°C)  at  S2  and  thick  black 
line  denotes  mixed  layer  depth.  Low-pass-filtered  (b)  east-west  U  and  (c)  north-south  V 
velocity  components  (cms-1)  at  M5  where  water  depth  is  47  m.  Cutoff  frequency  of  low- 
frequency  currents  is  0.08  cph.  Black  lines  indicate  MLD  at  S2. 


2011.  The  MLD  was  estimated  from  moored  tempera¬ 
ture  time  series  at  S2-S4,  because  salinity  had  a  negligi¬ 
ble  impact  here  on  the  vertical  structure  of  density  (TG). 
The  MLD  was  computed  as  the  depth  at  which  tem¬ 
perature  has  decreased  by  0.2°C  from  the  shallowest 
temperature  measurement.  Estimates  of  MLD  at  M2- 
M4  are  similar,  and  therefore  the  MLD  at  M2  is  used  for 
this  analysis  (Fig.  4a).  Note  that  there  was  not  a  TCP 
string  mooring  at  M5,  and  therefore  MLD  is  inferred 
from  the  temperature  records  at  M2-M4  (Fig.  4a).  After 
28  October  (year  day  300),  the  mixed  layer  around  the 
bank  deepened  below  47  m,  indicating  that  the  water 
column  over  the  bank  was  well  mixed.  A  detailed  dis¬ 
cussion  of  MLD  and  its  relationships  with  longer  time 
scale  processes  are  described  by  TG.  In  the  present 
study,  we  use  MLD  as  a  scaling  parameter  that  controls 
the  vertical  distribution  of  high-frequency  velocity 
fluctuations  in  the  mixed  layer. 

TG  described  low-frequency  flows  around  the  bank  for 
2011.  They  note  that  background  low-frequency  motions 
were  dominated  by  semidiurnal  tides  (12.42  h),  near- 
inertial  waves  (25.58  h),  and  60-  and  96-h  oscillations.  In 
general,  barotropic  diurnal  and  semidiurnal  tides  were 
weak,  and  mean  surface  currents  were  eastward  in  the 
upper  80  m.  Mean  currents  below  approximately  20-m 
above  the  bottom  followed  the  bathymetry.  Figures  4b 


and  4c  show  low-frequency  horizontal  currents  at  M5. 
Near-inertial  currents  were  typically  stronger  and  highly 
baroclinic  during  summertime  when  stratification  was 
strong.  Once  the  stratification  was  broken  down  during 
fall  and  winter  seasons,  the  currents  became  nearly  baro¬ 
tropic  with  relatively  weak  shear  in  the  water  column. 

4.  High-frequency  velocity  fluctuations 

Vertical  velocity  w  in  the  mixed  layer  over  the  bank 
was  as  large  as  10  cm  s-1  during  wind  events  (Fig.  2e);  w 
is  based  on  15-min-averaged  ADCP  measurements. 
Vertical  velocity  became  strongest  when  wind  speed  and 
Stokes  drift  velocity  were  largest  (Figs.  2a-d).  The 
squared  vertical  velocity  averaged  over  water  depth  w2H 
clearly  correlated  with  the  occurrence  of  surface  waves 
over  the  bank.  Figures  2d  and  2e  further  illustrate  the 
near-perfect  coincidence  between  the  vertical  velocity 
at  5.5  m  and  the  surface  wave  signature,  measured  by  the 
bottom-mounted  Ppod.  A  total  of  12  wind  events  with 
significantly  large  w2H  were  identified  between  June  and 
December  2011  (Fig.  2f).  In  general,  the  production- 
scale  turbulent  flows  in  the  mixed  layer  have  time  scales 
varying  from  several  minutes  to  hours.  These  time  scales 
are  determined  by  the  depth  of  the  mixed  layer  and  the 
dominant  turbulent  eddies  in  the  mixed  layer.  Therefore, 
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Fig.  5.  Variance-preserving  frequency  spectra:  (a)  u,  (b)  v,  and  (c)  vertical  w  velocity  components  computed  at  5.5-,  11. 5-,  17.5-,  and  23.5-m 

water  depths  for  the  wind  event  4.  The  inertial  frequency  is  0.041  cph. 


the  15-min  sampling  rate  may  partially  resolve  production- 
scale  turbulent  flows  in  the  mixed  layer.  We  suspect  that 
large  vertical  velocities  over  the  bank  were  related  to 
wave-induced  currents,  such  as  Langmuir  flows.  Ac¬ 
cording  to  the  theory  of  Craik  and  Leibovich  (1976),  the 
combined  effects  of  Stokes  drift  and  background  vor- 
ticity  IX  of  mean  currents  generates  a  vortex  force  F, 
which  drives  the  Langmuir  circulation,  where 

F  =  Us(z)xfl.  (3) 

In  the  following,  we  examine  these  high-frequency  mo¬ 
tions  in  the  mixed  layer  and  their  relationships  to  winds 
and  surface  waves. 

a.  Three-dimensional  velocity  structure 

To  examine  strength  and  distribution  of  velocity  vari¬ 
ance  at  higher  frequencies,  velocity  spectra  were  com¬ 
puted  during  several  wind  events.  One  example,  shown  in 
Fig.  5,  is  the  variance-preserving  velocity  spectra  for 
event  4.  A  spectrum  for  a  given  depth  was  computed  from 
512  data  points  (or  5.33-day  record)  centered  on  the  wind 
event.  The  spectra  were  multiplied  by  the  corresponding 
frequencies  to  construct  the  variance -preserving  spec¬ 
tra,  which  were  then  averaged  into  32  frequency  bands 
(Fig.  5).  Depth-time  series  of  horizontal  and  vertical 
velocities  for  event  4  are  plotted  in  Fig.  6.  Low-frequency 
currents  and  high-frequency  velocity  variances  are  illus¬ 
trated  separately.  Low-frequency  currents  were  con¬ 
structed  by  removing  velocities  at  higher  than  0.08  cph, 
while  high-frequency  velocity  currents  were  constructed 
by  removing  velocities  at  frequencies  lower  than  0.2  cph. 
A  fourth-order  Butterworth  filter  was  used  for  filtering 


the  velocity  fields.  Velocity  spectra  show  distinct  features: 
(i)  spectral  peak  in  horizontal  velocity  spectra  at  near- 
inertial  frequency  (0.041  cph),  and  (ii)  elevated  velocity 
variance  at  higher  frequencies  between  0.2  cph  and  the 
Nyquist  frequency  of  2  cph  (Fig.  5).  The  high-frequency 
energy  levels  were  largest  near  the  Nyquist  frequency, 
indicating  that  the  15-min  sampling  rate  would  resolve 
higher-frequency  motions  partially. 

During  event  4,  winds  increased  from  4  to  15 ms-1; 
the  southerly  winds  rotated  counter  clockwise  to  form 
westerly  winds  in  the  middle  part  of  event  (days  246- 
247)  before  veering  toward  the  south  during  the  latter 
part  of  the  event  (Fig.  2).  The  mixed  layer  deepened 
from  20  to  40  m;  and  low-frequency  currents  (such  as 
near-inertial  waves)  as  large  as  0.5 ms-1  were  formed. 
Here,  (u2)  was  smaller  than  (u2),  but  the  strongest  values 
of  (u2)  and  (v2)  were  in  the  crosswind  direction  (Figs.  2 
and  6);  (w2)  varied  with  winds,  but  followed  Stokes  drift 
closely  (not  shown).  The  high-frequency  velocity  vari¬ 
ances  ((u2),  (v2),  and  (w2))  decreased  rapidly  with  depth; 
for  example,  a  couple  of  orders  of  magnitude  difference  in 
variance  of  (w2)  were  found  between  5.5-  and  40-m  water 
depths  (Fig.  6).  Here,  (w2)  penetrated  to  greater  depths 
than  (u2)  and  (v2).  At  the  sea  surface  w  must  be  zero 
and  has  a  subsurface  maximum  between  the  surface 
and  5.5  m,  because  w  decays  with  depth  starting  at  5.5-m 
depth.  Furthermore,  high-frequency  motions  were  strong 
when  near-inertial  waves  were  large  (Fig.  6).  We  suspect 
that  the  advection  of  high-frequency  velocity  variances  by 
near-inertial  currents  might  be  an  important  term  in  the 
energy  balance  of  high-frequency  motions. 

We  further  examined  other  wind  events  to  understand 
the  three-dimensional  velocity  structure  in  the  mixed 
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Fig.  6.  Depth-time  sections  of  currents  at  M5  during  wind  event  4.  (left)  Low-pass  filtered  currents  (a)  U,  (c)  V, 
and  (e)  W  (ms-1).  East,  north,  and  upward  directions  are  positive,  (right)  Logio  values  of  squared  high-frequency 
(denoted  by  angle  brackets)  currents  (ms-1)  (b)  log10((z/2)),  (d)  log10((u2)),  and  (f)  log10((w2)).  The  black  lines 
denote  the  MLD.  The  cutoff  frequency  of  high-frequency  fluctuations  is  set  to  0.2  cph. 


layer.  Figure  7  illustrates  two  examples  of  depth-time 
variability  of  ( u 2),  ( v 2),  and  (w2),  and  magnitudes  and 
directions  of  Stokes  drift  and  surface  currents.  During 
the  summer  event  (event  2)  the  winds/Stokes  drifts  were 
approximately  from  the  south  and  the  mixed  layer  was 
shallow  (—20  m)  (Figs.  7a-d);  between  year  day  178  and 
179,  mixed  layer  low-frequency  currents  were  approxi¬ 
mately  90°  out  of  phase  with  Stokes  drift;  after  year  day 
179,  mixed  layer  currents  were  dominated  by  clockwise- 
rotating,  near-inertial  motions,  and  increasing  mean 
flow  (Fig.  7e).  Stokes  drift  of  about  3-5  cm  s-1  was  from 
the  southwest,  and  reduced  to  below  2  cm  s-1  after  year 
day  181  (Fig.  7a).  At  the  beginning  of  event  2,  both 
the  north-south  (crosswind)  velocity  variance  (i.e.,  (u2)) 
and  the  east-west  (downwind)  velocity  variance  (i.e.,  (v2)) 
were  generated.  The  vertical  velocities  were  strong 
and  (iv2)  strengthened,  penetrated  below  the  MLD  as 
the  Stokes  drift  increased  to  above  a  threshold  value 
(— 2cms_1),  and  diminished  as  Stokes  drift  decreased  be¬ 
low  the  threshold  (Fig.  7e).  However,  there  is  no  apparent 
relationship  between  the  magnitude  of  high-frequency 


motions  and  the  direction  of  near-inertial  currents  (Figs. 
7c,e,g,i).  During  the  fall  event  (event  5),  the  winds/ 
Stokes  drifts  were  approximately  from  the  east  and 
the  mixed  layer  was  deep  (— 40  m).  The  direction  of  the 
surface  current  was  nearly  180°  out  of  phase  with  the 
Stokes  drift  between  year  day  278  and  281  and  those 
currents  became  in  phase  during  the  later  stages  (Fig. 
7d).  At  the  beginning  of  event  5,  winds/Stokes  drift  were 
small;  Stokes  drift  opposed  the  surface  current;  and 
high-frequency  velocity  components  were  small.  During 
the  latter  part  of  event  5,  higher  east-west  (downwind) 
(i.e.,  (u2))  and  vertical  velocity  variances  were  generated 
as  winds/Stokes  drift  increased  and  both  Stokes  drift  and 
the  surface  current  were  in  the  same  direction  (coming 
from  the  east)  (Fig.  7d),  but  the  north-south  (crosswind) 
surface  current  component  (i.e.,(u2))  remained  low 
(Fig.  7j).  For  a  typical  Langmuir  flow  in  the  open-ocean 
surface  layer,  we  expect  generation  of  crosswind  surface 
currents,  but  here  we  observed  a  strong  downwind  sur¬ 
face  current  (Fig.  7h),  perhaps  similar  to  shear-driven 
turbulence.  These  examples  show  the  generation  of 
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Year  Day  Year  Day 

Fig.  7.  (a),(b)  The  magnitude  of  Us0  (red)  and  surface  current  C/c/5  (blue)  in  the  upper  15.5  m  at  M5  (cms-1) 
during  (left)  event  2  and  (right)  event  5.  (c),(d)  Corresponding  directions  (°)  of  Stokes  drift  0us  (red)  and  surface 
current  6uc  (blue).  (e)-fl)  Corresponding  depth-time  sections  of  squared  high-frequency  currents  (w2),  (u2),  and  ( v 2). 
Logio  values  of  velocity  (ms-1)  variances  are  plotted.  Thick  black  lines  denote  MLD.  Directions  follow  convention 
used  in  meteorology  and  represent  the  direction  from  which  the  Stokes  drift  and  currents  originate. 


a  three-dimensional  complex  flow  field  in  the  presence 
of  wind  events  and  background  currents.  However,  the 
timing  and  the  magnitude  of  velocity  variances  appear 
to  depend  on  both  strengths  and  directions  of  winds  and 
waves  over  the  bank.  We  further  examine  high  wind 
events  and  their  parameter  dependence  on  (w2)  in 
section  4e. 

b.  Probability  distribution 

As  discussed  above,  high-frequency  flows  observed 
during  wind  events  were  strongest  near  the  surface  and 


decayed  rapidly  with  depth.  We  also  noted  that  magni¬ 
tudes  of  upward  and  downward  motions  were  compa¬ 
rable.  Skewness  of  n  for  the  12  events  between  5.5  and 
15.5  m  is  -0.34  with  a  standard  deviation  0.36,  and  the 
corresponding  kurtosis  (or  flatness)  is  4.67  with  a  stan¬ 
dard  deviation  0.23.  These  statistics  show  that  the 
probability  distribution  of  w  differs  from  the  Gaussian 
distribution  (in  which  skewness  is  0  and  kurtosis  is  3). 
Figure  8  shows  cumulative  probability  distributions  of 
normalized  u ,  u,  and  w  and  the  corresponding  cumula¬ 
tive  normal  Gaussian  distributions  at  5.5-m  depth  for 
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Fig.  8.  Cumulative  probability  distribution  function  Cpdf  of  high-frequency  velocity  at  5.5-m  depth  is  plotted 
against  the  normalized  velocity  cf)  for  (a)  summer  event  2,  (b)  fall-winter  event  12.  The  black  triangles,  blue  open 
circles,  and  red  dots  denote  u,  v,  and  w  velocity  components,  respectively.  The  red  dashed  line  represents  the 
Gaussian  distribution  function. 


two  wind  events  (2  and  12),  where  the  normalized 
variable  </>  =  (O  -  0)/cr$,  O  represents  u ,  v ,  and  w,  cr®  is 
the  standard  deviation,  and  the  overbar  denotes  the 
mean  of  the  data  record.  The  events  shown  in  Fig.  8 
characterize  two  different  background  stratifications 
and  wind/wave  patterns  (Figs.  2  and  4).  Cumulative 
distributions  of  u ,  u,  and  w  are  similar,  and  they  all 
deviate  from  the  Gaussian  distribution  when 
which  is  nearly  20%  of  the  data  record.  We  noted 
similar  velocity  distributions  at  different  depths  within 
the  mixed  layer.  In  a  turbulent  flow,  the  values  of  flat¬ 
ness  are  large  relative  to  a  Gaussian  distribution  if  the 
probability  distribution  function  has  relatively  large 
values  in  its  tails  (e.g.,  Fig.  8).  The  flatness  occurs  when 
the  time  series  of  the  velocity  field  contains  significant 
numbers  of  sharp  peaks.  We  suspect  that  these  high- 
frequency  peaks  represent  production-scale  wind- 
driven  turbulent  eddies  and  coherent  structures  such  as 
Langmuir  circulation  cells. 

c.  Spatial  variability 

Figure  9  illustrates  high-frequency  kinetic  energy  at 
M1-M5  for  three  wind  events  (2,  4,  and  12).  Vertical 
profiles  of  horizontal  kinetic  energy  (HKE)  and  vertical 
kinetic  energy  (VKE)  were  computed  from  high-pass- 
filtered  AD  CP  velocity  profiles.  Kinetic  energy  profiles 


for  these  selected  events  were  constructed  by  time  av¬ 
eraging  of  squared  high-pass-filtered  velocities  ( u ,  v ,  w), 
where  VKE  =  0.5«w2»  and  HKE  =  0.5(«w2»  +  «v2»)  and 
«»  denotes  a  time  average  over  wind  events.  Here,  the 
cutoff  frequency  was  set  to  0.2  cph.  As  we  mentioned 
earlier,  these  events  (Fig.  9)  represent  different  back¬ 
ground  stratification  and  wind/wave  patterns  (Figs.  2 
and  4).  In  event  2,  the  MLD  was  shallow  (—20  m)  and 
winds  were  from  the  south-southeast  at  an  average 
speed  of  about  8 ms-1.  During  event  4,  the  MLD 
deepened  from  20  to  40  m  and  winds  accelerated  from  4 
to  15 ms-1  while  wind  and  waves  rotated  clockwise 
from  north  to  south.  During  event  12,  there  was  a  frontal 
passage  and  the  water  column  was  well  mixed.  The  high- 
frequency  kinetic  energy  at  M5  was  the  largest  of  all 
mooring  observations.  The  differences  in  energy  were 
more  pronounced  in  the  vertical  kinetic  energy  profiles 
than  in  the  horizontal  kinetic  energy  profiles.  A  signifi¬ 
cantly  high  energy  level  at  M5  could  be  caused  by  the 
intensification  of  surface  waves  over  the  bank.  M5  is 
located  at  47-m  water  depth  and  3-km  north  of  the  peak, 
which  is  about  20-m  below  the  sea  surface.  The  region 
north  of  the  peak  is  a  plausible  wave-focusing  region  for 
waves  propagating  from  the  south.  We  will  discuss  these 
wave-focusing  effects  and  their  impacts  on  Stokes  drift 
and  vertical  velocity  in  section  5. 
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Fig.  9.  (a)-(c)  Vertical  profiles  of  VKE  and  (d)-(f)  HKE  averaged  over  three  wind  events  [(top)  2,  (middle)  4,  and 
(bottom)  12].  Thick  solid  lines  represent  energy  levels  at  M5  and  dashed  lines  show  95%  confidence  limits.  Energy 
levels  at  Ml,  M2,  M3,  and  M4  are  marked  in  crosses,  open  circles,  open  diamonds,  and  open  triangles,  respectively. 
Confidence  limits  of  energy  levels  at  M1-M4  are  not  plotted. 


d.  The  e-folding  scales 

The  entrainment  depth  of  vertical  velocity  can  be 
examined  from  the  e-folding  length  of  tv2,  because  the 
e-folding  length  compares  individual  profiles  indepen¬ 
dently  of  their  magnitudes.  The  e-folding  length  £  of  tv2 
was  estimated  by  fitting  the  exponential  profile  between 
the  shallowest  depth  (5.5  m)  and  30.5  m,  provided  that  w 
at  the  shallowest  depth  was  greater  than  3  cms-1.  This 
minimizes  uncertainties  of  £  by  maintaining  well- 
resolved  vertical  profiles  for  the  computation.  For  the 
Craik  and  Leibovich  (1976)  theory,  the  Langmuir  cir¬ 
culation  is  governed  by  the  vortex  force  [Eq.  (3)],  where 
the  magnitude  and  vertical  structure  of  w2  depend 
strongly  on  Stokes  drift  [i.e.,Uy(z)].  Stokes  drift  decays 
exponentially  with  depth  and  is  a  function  of  horizontal 
wavelength  determined  by  wave  period  and  water 
depth.  The  e-folding  length  17  of  |Us(z)|  was  also  de¬ 
termined  by  fitting  an  exponential  profile  between  the 


surface  and  15-m  water  depth,  provided  the  magnitude 
of  the  Stokes  drift  at  the  surface  (i.e.,  Us0)  is  greater  than 
lcms”1.  The  e-folding  length  of  |Us(z)|  was  also  esti¬ 
mated  by  finding  a  depth  at  which  |Us(z)|/f/s0  =  e-1 
(e.g.,  Harcourt  and  D’Asaro  2008).  Both  methods  pro¬ 
duced  comparable  estimates  of  e-folding  lengths.  Note 
that  w2  profiles  were  available  every  15  min  while  Stokes 
drift  profiles  were  available  every  two  hours. 

Figure  10  illustrates  time  series  of  e-folding  lengths  £ 
and  r]  along  with  the  mixed  layer  depth  at  S2.  The 
e-folding  length  of  w2  varied  between  3.5  and  45  m 
(Fig.  10b)  with  an  average  and  a  standard  deviation  of 
8.3  and  3.9  m,  respectively  for  the  12  events.  Large 
values  of  £  were  found  with  deep  MLD  (Figs.  10a, b). 
Because  the  vertical  scale  of  turbulent  eddies  depends 
strongly  on  the  depth  of  the  mixed  layer,  we  can  expect 
large  values  of  £  for  deep  mixed  layers.  The  e-folding 
length  of  U$(z)  varied  between  2  and  9  m  (Fig.  10c)  with 
mean  rj  of  4.45  and  standard  deviation  of  1.02  m,  which  is 


August  2013 


WI JESEKER A  ET  AL. 


1639 


J _ I _ I _ I _ I _ I _ I _ I _ I _ I _ 

160  180  200  220  240  260  280  300  320  340 


Year  Day 

Fig.  10.  (a)  Time  series  of  MLD  at  S2.  Black  dashed  line  denotes  water  depth  (47  m)  at  M5. 
(b)  The  e-folding  length  (i.e.,  £)  of  w2  at  M5  for  individual  profiles,  (c)  The  e-folding  length  of 
Stokes  drift  (i.e.,  17)  (black  lines)  and  the  vertical  scale  of  the  peak  wave  [1/(2 kp);  red  lines]  at 
M5.  Numbers  1-12  are  the  wind  events  marked  in  Fig.  2. 


much  less  than  the  mixed  layer  depth  (Fig.  10a). 
Therefore,  wave  forcing  by  the  vortex  force  [Eq.  (3)]  is 
limited  to  the  upper  10  m  or  less.  As  shown  in  Fig.  10c, 
where  kp  is  the  peak  wavenumber,  the  vertical  scale  of 
dominant  surface  wave  (2/q,)-1  is  comparable  with  77 
indicating  that  the  vertical  structure  of  the  Stokes  drift 
is  controlled  by  the  dominant  wave.  The  mean  and 
standard  deviation  of  (2/c^)-1  for  the  12  events  are  4.75 
and  1.37  m,  respectively.  The  normalized  77  by  the 
mean  of  (2/q,)-1  is  about  0.94,  which  is  approximately 
a  factor  of  5  larger  than  the  corresponding  ratio  for  the 
PM  spectral  model  (Li  and  Garrett  1993).  Unlike  £,  77 
did  not  change  with  mixed  layer  depth  but  varied  with 
wind  speed  (Figs.  2  and  10).  On  average,  the  vertical 
kinetic  energy  penetrated  approximately  twice  the 
average  penetration  depth  of  the  Stokes  drift.  Our 
observations  indicate  that  the  penetration  depth  of  w  is 
determined  by  several  factors,  including  Us0,  £/*, 
MLD,  and  mean  currents,  while  the  penetration  of  the 
vortex  force  is  determined  by  the  e-folding  length  of 
the  Stokes  drift. 


e.  Scaling  vertical  velocity 

The  relative  importance  of  wind  and  wave  forcing  is 
described  by  the  turbulent  Langmuir  number  (McWilliams 
et  al.  1997;  Li  et  al.  2005;  Harcourt  and  D’Asaro  2008), 

La  t  =  (UJUs)m.  (4) 

t  J0 

Figure  11  shows  Laf  versus  ( w2D)IU% ,  where  (w2D)  was 
constructed  by  averaging  w2  over  the  MLD  and  then 
averaging  temporally  over  a  2-h  segment  before  dividing 
by  U%.  Lar  is  based  on  two-hourly  estimates  of  Us0  and 
U *.  Almost  all  of  the  wind  events  had  Laf  <  1,  although 
(wp)/Ul  varied  between  0.1  and  60.  Our  estimate  of 
averaged  Laf  for  the  12  wind  events  is  about  0.63  with 
a  standard  deviation  of  0.17  (Fig.  11)  and  is  comparable 
to  the  estimates  of  Kukulka  et  al.  (2011)  based  on 
measured  directional  wave  spectra  with  fc  =  0.4  Hz.  For 
a  fully  developed  wave  field  represented  by  the  PM 
spectral  model,  Laf  is  about  0.3  (Li  and  Garrett  1993;  Li 
et  al.  2005).  The  turbulent  Langmuir  number  is  expected 
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Fig.  11.  (wp)/Ul  vs  turbulent  Lar.  Black  dots  represent  the  entire  dataset  and  circles  represent 
data  from  12  wind  events  shown  in  Fig.  2. 


to  be  higher  than  0.3  for  developing  seas  under  fetch- 
limited  conditions  as  a  result  of  weaker  Stokes  drift  (Li 
et  al.  2005;  Harcourt  and  D’Asaro  2008). 

Another  important  scaling  parameter  is  the  Hoenik- 
ker  number  (Li  et  al.  2005),  which  describes  the  relative 
importance  of  convection  (unstable  buoyancy  forcing) 
and  vortex  force  driving  Langmuir  circulation.  How¬ 
ever,  we  do  not  have  the  net  surface  buoyancy  flux  to 
evaluate  the  Hoenikker  number.  We  examined  day  and 
night  variability  of  (w2)  for  cases  where  La t  <  1  and 
(■ w2d)IU%  >  1,  but  did  not  find  significant  differences 
during  wind  events. 

The  vertical  velocity  is  examined  within  the  frame¬ 
work  of  the  wind/wave  mechanism  such  as  Langmuir 
circulation  (Craik  and  Leibovich  1976).  This  framework 
is  often  used  to  describe  the  surface  velocity  variance 
observed  from  Doppler  sonar  systems  (e.g.,  Plueddemann 
et  al.  1996;  Smith  1998)  and  the  large-eddy  simulation  of 
the  oceanic  mixed  layer,  where  the  vortex  forcing  is  added 
to  mimic  Langmuir  circulation  (Skyllingstad  and  Denbo 
1995;  McWilliams  et  al.  1997;  Harcourt  and  D’Asaro  2008; 
Grant  and  Belcher  2009).  D’Asaro  and  Dairiki  (1997) 
and  Tseng  and  D’Asaro  (2004)  compared  MLD  aver¬ 
aged  and  squared  vertical  velocity  with  squared  friction 
velocity,  and  found  that  on  average,  {w2D)  was  about 
a  factor  of  1. 5-3.0  larger  than  showing  the  departure 

from  the  wall  boundary  layer  theory.  Their  empirical 
evidence  is  somewhat  surprising  because  of  lack  of 


dependence  on  surface  wave  parameters.  Plueddemann 
et  al.  (1996)  suggest  that  surface  velocity  variance  scales 
as  U*Us0,  while  Smith  (1998)  reports  that  velocity  var¬ 
iance  follows  t/|o  scaling  once  Langmuir  circulation  is 
established.  Harcourt  and  D’Asaro  (2008)  suggest  that 
turbulent  vertical  velocity  scales  as  £/* La^2/3,  where 
La l  =  (U */ (Usl))112 ,  and  USl  is  the  averaged  Stokes 
drift  within  the  surface  layer  L,  which  is  taken  as  20%  of 
the  height  of  the  mixed  layer.  Based  on  LES  studies, 
Grant  and  Belcher  (2009)  report  that  the  velocity  scale 
of  Langmuir  turbulence  can  be  expressed  as  (U*Us0)113, 
similar  to  the  scaling  suggested  by  Smith  (1996).  The 
observations  and  modeling  studies  indicate  that  there  is 
no  consistent  scaling  for  Langmuir  cell  velocities,  even 
though  all  these  scaling  arguments  are  either  related  to 
Us0  or  a  combination  of  17* ,  Us0,  and  mixed  layer  depth. 
Here,  (w2)  can  vary  with  surface  forcing  such  as  Stokes 
drift,  friction  velocity,  and  surface  buoyancy  flux,  and 
with  upper-ocean  dynamics  such  as  mixed  layer  depth 
and  advection  caused  by  mean  currents. 

We  will  apply  some  of  the  existing  scaling  methods  to 
the  observed  vertical  velocity  in  the  mixed  layer  over  the 
bank.  We  limit  our  analysis  to  12  wind  events  shown  in 
Fig.  2f  where  a  significant  number  of  good  vertical  ve¬ 
locity  profiles  can  be  found.  These  wind  events  occupied 
a  total  of  about  41  days  between  June  and  December 
2011.  We  used  3928,  15-min-averaged  w  profiles  for 
the  analysis.  We  note  that  the  magnitude  and  vertical 
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Fig.  12.  {wp)/Ul  vs  fractional  MLD  D/H,  where  A,  B,  and  C 
represent  three  groups  of  (wp)/Ul,  and  each  group  is  divided  into 
four  subgroups  of  D/H  as  indicated  by  1-4.  The  colors  represent 
different  D/H  groups. 

structure  of  w2  and  the  magnitudes  and  directions  of 
Stokes  drifts  and  winds  varied  significantly  within  a  wind 
event.  Therefore,  the  scaling  of  vertical  profiles  should 
be  conducted  by  identifying  profiles  that  have  similar 
dynamical  regimes  or  belong  to  similar  parameter 
regimes  defined  by  the  underlying  physical  process. 
Figure  12  shows  spreading  of  (wfy/U*  as  a  function  of 
D/H,  where  D/H  is  the  fractional  mixed  layer  depth.  The 
segments  cover  a  broad  range  of  values  where  {w2D)IU% 
ranges  from  0.1  to  60  and  D/H  ranges  from  0.2  to  1 
(Fig.  12).  The  vertical  velocity  profiles  were  divided 
into  three  groups  of  (wfy/Ul  denoted  by  A,  B,  and  C 
(Fig.  12).  These  groups  were  selected  by  considering  the 
mixed  layer  vertical  kinetic  energy  relative  to  the  wall 
boundary  layer  estimates;  for  example,  (wfy/U*  >  1 
represents  regimes  which  depart  from  wall  boundary 
scaling.  As  shown  in  Fig.  12,  the  data  belong  to  a  wide 
range  of  mixed  layer  depths.  Therefore,  each  group  was 
divided  into  four  subgroups  representing  four  different 
segments  of  D/H  (e.g.,  Al,  A2,  A3,  and  A4  in  Fig.  12).  For 
example,  the  subgroup  Al  represents  weak  energy  and 
shallow  mixed  layer  profiles  whereas  the  group  A4  rep¬ 
resents  weak  energy  and  deep  mixed  layer  profiles.  Profiles 
of  (w2)  in  each  subgroup  were  scaled  by  U%,  U*Us0 , 
( UlUs0 )2,  and  Uj()  before  constructing  averaged  profiles. 

Scaled  (w2)  profiles  of  all  subgroups  are  compared  as 
a  function  of  scaled  depth,  z/D  (Fig.  13).  We  compare 
near-surface  estimates  and  the  vertical  structure  of 
scaled  (w2)  resulting  from  U%,  U*Us0 ,  and  £/|  ,  because 
the  scaling  by  U*Us0  and  ( U%UsQ )2  produce  relatively 


similar  results.  The  averaged  profiles  of  (w2),  scaled  by 
Ul,  U*Us0 ,  and  UjQ,  are  plotted  in  Figs.  13a-c,  Figs. 
13d-f,  and  Figs.  13g-i,  respectively.  These  contain  pro¬ 
files  from  groups  A,  B,  and  C  (see  Fig.  12),  respectively. 
There  is  a  fair  qualitative  agreement  with  U%  scaling  for 
group  A  (Fig.  13a).  Here,  (w2)/Ul  profiles  in  group  A 
collapse  to  a  single  curve  in  the  upper  half  of  MLD. 
However,  U%  scaling  fails  for  (w2)/Ul  >  1  (Figs.  13b, c). 
The  scaling  of  (w2)  by  U*Us0  is  also  qualitatively  con¬ 
sistent  with  group  A  even  though  near-surface  estimates 
are  slightly  smaller  than  the  expected  estimate  of  unity 
(Fig.  13d).  Scaling  of  U*  Us0  over  predicts  slightly  for  Bl, 
B3,  and  B4  (Fig.  13e).  Like  U%  scaling,  U*Us0  does  not 
provide  a  satisfactory  agreement  for  group  C,  where 
(■ w2)/U*Us0  near  the  surface  is  about  10.  There  is  a  wide 
range  of  values  for  (w2)/U£q  for  groups  A  and  B,  but 
(w2)/Ujo  near  the  surface  is  close  to  unity  for  subgroups 
A2,  A4,  Bl,  B2,  and  B4.  None  of  the  scaling  methods 
works  well  for  group  C,  nevertheless  scaling  pro¬ 
vides  near-surface  estimates  of  (w2)/Ujo  close  to  a  factor 
of  2.  In  most  of  the  groups,  the  vertical  structure  of 
scaled  profiles  decays  with  varying  e-folding  scales. 
These  results  reflect  the  fact  that  there  is  no  unique 
scaling  of  (w2)  profiles  in  the  mixed  layer  based  on 
surface  forcing  parameters,  although  (w2)  depends  on 
t/*,  Us0,  mixed  layer  depth,  and  e-folding  scales  of  w2. 

Wind  events  generate  near-inertial  currents  (e.g.,  Fig. 
5),  and  the  relative  direction  between  winds  and  surface 
currents  can  make  a  difference  in  the  vertical  shear  of 
horizontal  currents  which  in  turn  can  alter  turbulence 
production  in  the  mixed  layer.  This  scenario  could  be 
a  plausible  explanation  why  our  parameterizations  were 
significantly  different  from  previous  publications.  We 
examined  the  impact  of  low-frequency  (near  inertial) 
currents  on  high-frequency  vertical  velocity  by  exam¬ 
ining  correlations  between  vertical  velocity  and  di¬ 
rectional  difference  between  winds  and  near-surface 
currents.  Relations  among  wind  and  surface  current  di¬ 
rections,  wind  speed,  and  variability  of  high-frequency 
flows  for  12  events  (Fig.  2)  were  examined.  We  did  not 
find  a  clear  relation  between  the  angular  difference  of  the 
mixed  layer  current  and  wind  speed,  and  (w2)  during 
different  stages  of  wind  events.  The  difference  between 
wind  and  current  directions  were  computed  from  2-h- 
averaged  current  and  wind  fields.  Here,  near-surface 
currents  at  M5  were  constructed  by  depth- averaging  low- 
pass-filtered  horizontal  currents  between  5.5  and  15.5  m. 

5.  Discussion 

a.  Vertical  scaling 

Scaling  of  (w2)  by  surface-forcing  parameters  such  as 
friction  velocity  and  Stokes  drift  produced  inconclusive 
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Fig.  13.  Scaled  squared  vertical  velocity  plotted  as  a  function  of  scaled  water  depth  (. zID )  (a)-(c)  U%  scaling  for  groups  A,  B,  and  C, 
respectively,  (d)-(f)  U *  Us0  scaling,  and  (g)-(i)  scaling  are  represented.  Colored  symbols  represent  regions  of  DIH  as  marked  in  Fig.  12. 


results.  We  note  that  the  e-folding  scale  of  (w2)  has 
a  wider  range  and  varies  with  the  mixed  layer  depth.  For 
example,  profiles  in  groups  B  and  C  have  higher  energy 
levels  at  depth  than  the  rest  of  the  profiles.  The  vertical 
distribution  of  (w2)  should  also  be  considered  when 
scaling  (w2)  (e.g.,  D’Asaro  and  Dairiki  1997;  Harcourt 
and  D’Asaro  2008),  because  individual  profiles  of  w2 
have  diverse  e-folding  scales  and  are  associated  with 
different  mixed  layer  depths  (Fig.  9).  Therefore,  we 
consider  mixed  layer-averaged  statistics  such  as  squared 
vertical  velocity  {w2D)  as  our  variable  to  be  scaled  and 
then  examine  the  variability  of  scaled  (w2D)  as  a  function 
of  the  e-folding  length  scale  of  w2  and  the  mixed  layer 


depth.  Figure  14  illustrates  (wfy/Ul  versus  f/Z)  and 
(wd)^^sq  versus  where  f/Z)  is  the  fractional  e-folding 
length.  For  (wfy/Ul  <  1,  the  scaling  by  U%  provides 
a  tighter  relationship  with  f/Z)  (Fig.  14a),  which  is  con¬ 
sistent  with  Fig.  13a,  but  U%  scaling  fails  for  {wp)IU%  >  1. 
The  scaling  by  UjQ  forces  close  fitting  between  (wp)/UjQ 
and  f/Z>.  As  shown  in  Fig.  14b,  log((w|,)/£/|0)  is  approx¬ 
imately  a  linear  function  of  f/Z),  and  the  least  squares  fit 
to  the  data  is 

^aExp^I) 


(5) 
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Fig.  14.  (a)  ( w2D)HJ%  vs  £/D  (normalized  e-folding  length  of  (w2)  by  mixed  layer  depth), 
(b)  (wp)/Uso  vs  £ID.  The  colored  symbols  represent  groups  of  fractional  mixed  layer  depth  DIH 
as  marked  in  Fig.  12.  The  thick  solid  line  in  (b)  represents  the  least  squares  exponential  fit. 


where  a  =  0.2,  /3  =  3.75,  and  0  <  fAD  <  1.  In  general,  UjQ 
scaling  provides  a  better  representation  of  data,  when 
effects  of  the  mixed  layer  depth  and  the  penetration 
depth  of  the  vertical  velocity  are  included. 

b.  Wave  focusing  over  the  bank 

We  noted  that  high-frequency  kinetic  energy  over  the 
bank  at  M5  is  significantly  larger  than  the  kinetic  energy 
around  the  bank,  and  these  differences  can  be  clearly 
seen  in  the  vertical  kinetic  energy,  especially  when  winds 
were  from  the  south-southeast  (Fig.  8).  We  suspect  that 
these  higher  velocity  fluctuations  are  generated  by 
wave-current  interaction  mechanisms  defined  by  the 
vortex  force  [Eq.  (3)].  The  amplitude  of  the  vertical 
component  of  the  vortex  force  is  Fv  =  \Us(z)dUm/dz\, 
where  |Uy(z)|  is  the  amplitude  of  the  Stokes  drift  and 
d  Umldz  is  the  vertical  gradient  of  the  mean  current  in  the 
direction  of  the  Stokes  drift.  To  examine  the  influence  of 
dUm/dz  on  Fv ,  these  gradients  are  compared  at  different 
sites.  As  shown  in  Fig.  15,  all  five  moorings  sites  over  and 
around  the  bank  had  similar  profiles  of  vertical  shear 
of  the  mean  currents.  Components  of  vertical  shear 
were  computed  from  centered  differencing  of  low-pass- 
filtered  ADCP  records,  where,  the  squared  shear  of  the 
mean  current  is  Sh2  =  ( dUm/dz )2  =  [( dU/dz )2  +  (dV/dz)2] 
(Fig.  15).  The  low-frequency  shear  is  dominated  by  near- 
inertial  currents  and  is  largest  at  the  base  of  the  mixed 


layer.  Shear  profiles  vary  with  wind  events  and  back¬ 
ground  stratification,  but  the  magnitude  and  vertical 
structure  of  shear  profiles  are  relatively  similar  for  all 
mooring  sites.  Therefore,  it  is  unlikely  that  dUm/dz  alone 
generates  a  larger  vertical  velocity  over  the  bank  than  at 
other  locations  around  the  bank. 

The  propagation  of  surface  waves  over  complex 
bottom  bathymetric  features  may  be  associated  with 
many  processes  including  shoaling,  refraction,  energy 
dissipation,  and  diffraction.  We  suspect  that  waves 
coming  from  the  south  (from  southwest  to  southeast) 
generate  a  focusing  region  behind  the  peak  of  the 
bank,  while  enhancing  surface  wave  amplitudes  and 
Stokes  drift.  To  examine  the  plausibility  of  wave  fo¬ 
cusing,  we  apply  a  refraction  and  diffraction  numer¬ 
ical  wave  model  (REF/DIF1)  developed  by  Kirby 
and  Dalrymple  (1994).  REF/DIF1  is  a  phase-resolving 
parabolic  refraction-diffraction  model  for  ocean  sur¬ 
face  wave  propagation. 

We  apply  REF/DIF1  to  a  domain  of  14  km  (north- 
south)  by  10  km  (east-west)  encompassing  the  EFGB 
(Fig.  1)  with  the  x  axis  as  the  offshore  boundary  for  in¬ 
coming  waves.  The  surface  Stokes  drift  for  a  mono¬ 
chromatic  wave  can  be  estimated  as  Us0  =  a2kco ,  where 
a  is  the  amplitude  of  the  monochromatic  wave,  k  is  the 
wavenumber,  and  co  is  the  radian  frequency.  Note  the 
wave  height  is  twice  the  wave  amplitude. 
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Fig.  15.  Squared  vertical  shear  (i.e.,  Sh2)  of  low-frequency  horizontal  currents  at  M1-M5  for  three  different  wind  events:  (a)  2,  (b)  4,  and 
(c)  12.  Bullets,  open  circles,  diamonds,  triangles,  and  thick  solid  line  denote  Ml,  M2,  M3,  M4,  and  M5,  respectively.  Thick  dashed  lines 
denote  mean  MLD  at  M2  mooring  and  thin  dashed  lines  denote  one  std  dev  from  the  mean. 


Figure  16  shows  simulation  results  of  wave  height  and 
Stokes  drift  for  a  representative  monochromatic  wave  of 
1-m  amplitude  (or  2-m  wave  height)  and  8-s  wave  pe¬ 
riod,  propagating  from  the  southern  boundary  into  the 
model  domain.  The  selected  wave  statistics  are  similar  to 
the  waves  statistics  observed  during  event  2  (Fig.  2).  The 
model  results  display  the  spatial  variation  of  wave  height 
and  Stokes  drift  induced  by  the  depth  changes  from  the 
presence  of  the  bank.  The  northward-propagating  wave 
focuses  over  a  narrow  area  behind  the  peak  of  the  bank 
while  generating  a  much  higher  wave  height  and  Stokes 
drift.  The  wave  height  in  the  focusing  area  reaches  3  m 
which  is  a  50%  increase  over  the  incoming  wave.  The 
surface  Stokes  drift  increases  to  12cms_1,  which  is  ap¬ 
proximately  a  factor  of  2.5  higher  than  that  of  the  in¬ 
coming  wave.  Several  sensitivity  REF/DIF1  simulations 


were  also  conducted  by  changing  the  incoming  wave 
direction  by  ±15°.  Similar  narrow  band  wave-focusing 
zones  were  found.  Out  of  five  moorings,  only  one  moor¬ 
ing  (M5)  is  in  the  vicinity  of  the  focusing  zone.  Our  sim¬ 
ulation  is  limited  to  a  monochromatic  wave  train  with 
a  single  direction  and  no  background  currents.  Never¬ 
theless  this  simulation  demonstrates  the  underlining 
physical  process.  It  is  likely  that  the  observations  of  high 
energy  levels  of  kinetic  energy  at  M5  would  be  a  result  of 
enhancement  of  Stokes  drift  caused  by  wave  focusing 
induced  by  the  bathymetric  variation  of  the  EFGB. 

c.  Impact  of  vertical  motion  on  corals 

The  EFGB  marine  sanctuary  is  a  biologically  diverse 
environment  and  supports  a  large  number  of  species 
including  corals.  Coral  reefs  are  complex,  dynamic,  and 
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Fig.  16.  Color  shading  shows  spatial  variations  of  (left)  wave  height  and  (right)  Stokes  drift  from  REF/DIF1 
simulations.  Thin  dashed  black  contours  show  the  bathymetry  of  EFGB.  Black  squares  are  the  locations  of  Barny 
moorings  (M1-M5)  and  M5  is  located  over  the  bank  (see  also  Fig.  1). 


sensitive  ecosystems.  Therefore  understanding  of  the 
hydrodynamics,  turbulence  characteristics,  mixing  of 
water  masses  over  banks,  and  transport  of  suspended 
sediment  are  important  for  the  coral  reef  system.  Wave- 
current  interactions  and  generation  of  deeply  pene¬ 
trating  vertical  motions  can  play  an  important  role  in 
mixing,  rapid  deepening  of  the  mixed  layer,  and  verti¬ 
cal  transports  of  near-surface  water  toward  coral  hab¬ 
itats  on  the  bottom.  Vertical  velocities  intensify  during 
wind  events,  while  wave  focusing  can  further  enhance 
wave-driven  Stokes  drift  thus  providing  an  efficient 
mixing  and  transport  mechanism  over  the  bank. 

6.  Conclusions 

The  main  objectives  of  this  study  are  as  follows:  (i)  to 
examine  surface  wave  effects  on  high-frequency  flows 
over  the  EFGB,  (ii)  to  quantify  spatial  and  temporal 
characteristic  of  these  high-frequency  flows,  and  (iii)  to 
examine  the  scaling  laws  of  vertical  velocity  based  on  the 
Langmuir  circulation  theory.  The  major  observational 
findings  and  inferences  of  this  study  are  the  following. 

•  High-frequency,  three-dimensional  flows  with  frequen¬ 
cies  ranging  from  0.2  to  2  cph  (Nyquist  frequency)  were 
found  in  the  mixed  layer  when  wind  speed  and  Stokes 
drift  were  large.  It  is  suggested  that  Stokes  drift  and 
mean  flow  interactions  such  as  Langmuir  circulation 
are  responsible  for  generating  the  observed  velocity 
fluctuations. 

•  Stokes  drift  over  the  bank  was  estimated  from  two- 
dimensional  directional  wave  spectra  from  a  600-kHz 


ADCP.  The  estimated  Stokes  drift  at  the  surface  (i.e., 
Us0)  varied  with  the  wind  speed  and  is  about  0.6% 
of  the  wind  speed  for  winds  greater  than  10ms_1. 
However,  our  estimates  are  a  factor  of  2  smaller  than 
the  estimates  based  on  the  PM  spectral  model. 

•  Vertical  velocity  fluctuations  showed  upward  and 
downward  motions  with  similar  magnitudes.  The  prob¬ 
ability  distribution  of  velocity  fluctuations  ( u ,  v ,  w) 
departs  from  the  Gaussian  distribution  especially  at  the 
tail  ends  of  the  distribution. 

•  Squared  vertical  velocity  w2  was  strongest  near  the 
surface  and  decayed  exponentially  with  depth,  and 
the  associated  e-folding  length  is  a  factor  of  2  larger 
than  the  e-folding  length  associated  with  Stokes 
drift. 

•  The  2-h  and  mixed  layer-averaged  (w2D),  normalized 
by  the  squared  friction  velocity  XJ%,  varied  from  0.1  to 
60.  The  (wp)/Ul  was  strongest  when  the  turbulent 
Langmuir  number  was  less  than  unity  and  the  mixed 
layer  was  shallow. 

•  Scaling  of  squared  vertical  velocity  profiles  was  exam¬ 
ined  as  functions  of  t/*,  Us0,  and  UjQ.  Although  w2  is 
likely  to  be  a  function  of  Stokes  drift,  t/*,  and  mixed 
layer  depth,  a  unique  scaling  for  profiles  of  w2  was  not 
found.  For  (wfy/Ul  <  5,  there  is  a  qualitative  agree¬ 
ment  with  different  scaling  for  wide  ranges  of  the 
fractional  MLD.  For  {wjf/Ul  >  5,  all  the  scaling 
techniques  underpredict  near-surface  estimates  of  w2. 
However,  when  we  consider  mixed  layer-averaged 
statistics  such  as  squared  vertical  velocity  (w2D)  as  the 
scaling  variable,  the  scaling  by  t/|  forces  tighter  re¬ 
lationships  between  (w^yt/J  and  f/D,  and  (w^)/I/|o 
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is  approximately  an  exponential  function  of  £AD.  In 
general,  £/|o  scaling  provides  a  better  representation  of 
(wp)  data,  when  effects  of  the  mixed  layer  depth  and 
the  penetration  depth  of  the  vertical  velocity  are 
included. 

•  Largest  velocity  variances  were  found  over  the  bank  at 
M5,  located  nearly  3  km  north  of  the  peak  of  the  bank. 
It  is  suggested  based  on  a  wave  refraction-diffraction 
model  (REF/D IF1)  that  the  increase  in  kinetic  energy 
over  the  bank  is  due  to  enhancement  of  Stokes  drift  by 
topography-induced  waves  focusing. 

•  The  lack  of  agreement  with  open-ocean  parameteri- 
zations  of  vertical  kinetic  energy  over  the  bank  is 
likely  due  to  the  enhancement  of  kinetic  energy 
resulting  from  localized  focusing  of  surface  waves. 

Wave-driven  velocity  fluctuations  play  an  important 
role  in  mixing  and  vertical  transports  and,  in  addition, 
can  strongly  impact  coral  habitats  over  the  EFGB.  Our 
study  reveals  that  such  velocity  fluctuations  over  the 
bank  are  closely  related  to  surface  wave  conditions  that 
can  have  strong  localized  variability  as  a  result  of  wave 
focusing  over  the  bank.  Our  findings  indicate  that  the 
scaling  of  wave-induced  mixed  layer  currents  over 
a  shelf  edge  bank  differs  from  the  open  ocean.  Recently, 
Belcher  et  al.  (2012)  examined  the  global  perspective  of 
Langmuir  turbulence  in  the  ocean  surface  boundary 
layer  using  LES  and  diagnostic  analysis  techniques. 
However,  the  lack  of  understanding  of  detailed  physical 
process  suggests  that  further  studies  are  needed  for  the 
generalization  of  wave-driven  mixing  processes  in  the 
coastal  environment. 
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